Ideal GLM-MHD - a new mathematical model for simulating astrophysical plasmas

Magnetic fields are ubiquitous in space. As there is strong evidence that magnetic fields play an important role in a variety of astrophysical processes, they should not be neglected recklessly. However, analytic models in astrophysical either do often not take magnetic fields into account or can do this after limiting simplifications reducing their overall predictive power. Therefore, computational astrophysics has evolved as a modern field of research using sophisticated computer simulations to gain insight into physical processes. The ideal MHD equations, which are the most often used basis for simulating magnetized plasmas, have two critical drawbacks: Firstly, they do not limit the growth of numerically caused magnetic monopoles, and, secondly, most numerical schemes built from the ideal MHD equations are not conformable with thermodynamics. In my work, at the interplay of math and physics, I developed and presented the first thermodynamically consistent model with effective inbuilt divergence cleaning. My new Galilean-invariant model is suitable for simulating magnetized plasmas under extreme conditions as those typically encountered in astrophysical scenarios. The new model is called the "ideal GLM-MHD" equations and supports nine wave solutions. The accuracy and robustness of my numerical implementation are demonstrated with a number of tests, including comparisons to other schemes available within in the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code FLASH. A possible astrophysical application scenario is discussed in detail

Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics

This article serves as a summary outlining the mathematical entropy analysis of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD equations as they are particularly useful for mathematically modeling a wide variety of magnetized fluids. In order to be self-contained we first motivate the physical properties of a magnetic fluid and how it should behave under the laws of thermodynamics. Next, we introduce a mathematical model built from hyperbolic partial differential equations (PDEs) that translate physical laws into mathematical equations. After an overview of the continuous analysis, we thoroughly describe the derivation of a numerical approximation of the ideal MHD system that remains consistent to the continuous thermodynamic principles. The derivation of the method and the theorems contained within serve as the bulk of the review article. We demonstrate that the derived numerical approximation retains the correct entropic properties of the continuous model and show its applicability to a variety of standard numerical test cases for MHD schemes. We close with our conclusions and a brief discussion on future work in the area of entropy consistent numerical methods and the modeling of plasmas

Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics

This article serves as a summary outlining the mathematical entropy analysis of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD equations as they are particularly useful for mathematically modeling a wide variety of magnetized fluids. In order to be self-contained we first motivate the physical properties of a magnetic fluid and how it should behave under the laws of thermodynamics. Next, we introduce a mathematical model built from hyperbolic partial differential equations (PDEs) that translate physical laws into mathematical equations. After an overview of the continuous analysis, we thoroughly describe the derivation of a numerical approximation of the ideal MHD system that remains consistent to the continuous thermodynamic principles. The derivation of the method and the theorems contained within serve as the bulk of the review article. We demonstrate that the derived numerical approximation retains the correct entropic properties of the continuous model and show its applicability to a variety of standard numerical test cases for MHD schemes. We close with our conclusions and a brief discussion on future work in the area of entropy consistent numerical methods and the modeling of plasmas

An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and Numerical Verification

The first paper of this series presents a discretely entropy stable discontinuous Galerkin (DG) method for the resistive magnetohydrodynamics (MHD) equations on three-dimensional curvilinear unstructured hexahedral meshes. Compared to other fluid dynamics systems such as the shallow water equations or the compressible Navier-Stokes equations, the resistive MHD equations need special considerations because of the divergence-free constraint on the magnetic field. For instance, it is well known that for the symmetrization of the ideal MHD system as well as the continuous entropy analysis a non-conservative term proportional to the divergence of the magnetic field, typically referred to as the Powell term, must be included. As a consequence, the mimicry of the continuous entropy analysis in the discrete sense demands a suitable DG approximation of the non-conservative terms in addition to the ideal MHD terms. This paper focuses on the resistive MHD equations: Our first contribution is a proof that the resistive terms are symmetric and positive-definite when formulated in entropy space as gradients of the entropy variables, which enables us to show that the entropy inequality holds for the resistive MHD equations. This continuous analysis is the key for our DG discretization and guides the path for the construction of an approximation that discretely mimics the entropy inequality, typically termed entropy stability. Our second contribution is a detailed derivation and analysis of the discretization on three-dimensional curvilinear meshes. The discrete analysis relies on the summation-by-parts property, which is satisfied by the DG spectral element method (DGSEM) with Legendre-Gauss-Lobatto (LGL) nodes. Although the divergence- free constraint is included in the non-conservative terms, the resulting method has no particular treatment of the magnetic field divergence errors, which might pollute the solution quality. Our final contribution is the extension of the standard resistive MHD equations and our DG approximation with a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM). As a conclusion to the first part of this series, we provide detailed numerical validations of our DGSEM method that underline our theoretical derivations. In addition, we show a numerical example where the entropy stable DGSEM demonstrates increased robustness compared to the standard DGSEM

Supporting strategic crew management at passenger railways—model, method and system

This paper reports the results of a joint project with a large railway company in Germany to build a decision support system for analyzing the consequences of timetable changes, modifications of break and working time regulations as well as changes in the cost structure on future crew needs. For that purpose we have developed a mathematical model of the underlying crew scheduling problem that respects all the organizational and technical constraints as well as labor regulations. We have implemented a Branch&Price based optimization system that is used to perform scenario analyses of future crew needs using medium-term timetable drafts as input data

A uniquely defined entropy stable matrix dissipation operator for high Mach number ideal MHD and compressible Euler simulations

We describe a unique averaging procedure to design an entropy stable dissipation operator for the ideal magnetohydrodynamic (MHD) and compressible Euler equations. Often in the derivation of an entropy conservative numerical flux function much care is taken in the design and averaging of the entropy conservative numerical flux. We demonstrate in this work that if the discrete dissipation operator is not carefully chosen as well it can have deleterious effects on the numerical approximation. This is particularly true for very strong shocks or high Mach number flows present, for example, in astrophysical simulations. We present the underlying technique of how to construct a unique averaging technique for the discrete dissipation operator. We also demonstrate numerically the increased robustness of the approximation

Evolution of fluid flow and carbonate recrystallization rates in deep-sea sediments of the Equatorial Pacific

Fluid flow and carbonate recrystallization rates of deep-sea sediments from eight locations in the Equatorial Eastern Pacific were determined by using delta Ca-44/40 values of pore water and corresponding sediments. The studied drill sites of IODP Exp. 320/321 are located along a transect of decreasing crustal age and reveal different characteristic pore water depth profiles. The younger sites show an overall isotopic equilibration with the sediment in the upper part of the sedimentary column. In the lower part, the delta Ca-44/40 of the pore water increases back to seawater-like values at the sediment/basalt interface, forming a bulge-shaped pore water profile. The magnitude of the delta Ca-44/40 pore water bulge decreases with increasing age of the oceanic crust and sediment cover, resulting in seawater-like delta Ca-44/40 values throughout the sedimentary column in the oldest Sites U1331 and U1332. These findings indicate a seawater-like fluid input from the underlying crust into the sediment. Thus, after sedimentation, carbonate recrystallization processes start to enrich the pore water in Ca-40, and after a time of carbonate recrystallization and cooling of oceanic crust, a flow of seawater-like fluid starts to move upwards through the sedimentary column, enriching the pore water with Ca-44. We established a carbonate recrystallization and fluid flow model to quantify these processes. Our determined carbonate recrystallization rates between 0.000013e((-t/15.5)) and 0.00038e((-t/100.5)) and fluid flow rates in the range of 0.42-19 m*Myr(-1) indicate that the fluid flow within the investigated sites of IODP Exp. 320/321 depends on the sedimentary composition and location of the specific site, especially the proximity to a recharge or discharge site of a hydrothermal convection cell

A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD

Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation of the discrete entropy. Present works in the field of entropy stable numerical schemes are concerned with thorough derivations of entropy conservative fluxes for ideal MHD. However, as we show in this work, if the dissipation operator is not constructed in a very specific way, it cannot lead to a generally stable numerical scheme. The two main findings presented in this paper are that the entropy conserving flux of Ismail & Roe can easily break down for certain initial conditions commonly found in astrophysical simulations, and that special care must be taken in the derivation of a discrete dissipation matrix for an entropy stable numerical scheme to be robust. We present a convenient novel averaging procedure to evaluate the entropy Jacobians of the ideal MHD and the compressible Euler equations that yields a discretization with favorable robustness properties. (C) 2016 Elsevier Inc. All rights reserved